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' Abstract 

o ■ 

OO ■ We show that if the solar and atmospheric data are both described by maximal 

, vacuum oscillations at the relevant mass scales then there exists a unique mixing matrix 

I for three neutrino flavors. The solution necessarily conserves CP and automatically 

OO ' implies that there is no disappearance of atmospheric I'e, consistent with indications 

. from the Super-Kamiokande experiment. We also investigate the consequences for 

I three-neutrino mixing if the solar and atmospheric oscillations exhibit mixing that is 



large but not maximal. For non-maximal mixing fg ^ Ut and <-> oscillations are 
predicted that may be observable in future long-baseline experiments. 
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Introduction. It is now understood that neutrino oscillations can describe the data on the 
solar neutrino deficit H, ^, the atmospheric neutrino anomaly ^, P], and the results from 



the LSND experiment 0, only if a sterile neutrino is introduced p|, 0, |10|, ITl|, |T2[. However 
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a conservative 

approach is to assume that oscillations need only account for the solar and atmospheric 
data; also, it is possible to introduce new lepton-fiavor changing operators with coefficients 
small enough to evade present exclusion limits, but large enough to explain the small LSND 
amplitude fl^. Since a comparison of the theoretical expectations |]15| with the Super- 
Kamiokande (SuperK) atmospheric data suggests maximal oscillations with mass-squared 
difference 5m\^^ ^ 5 x 10~^ eV^ |p, and in light of the recent data from Super-Kamiokande 
favoring vacuum long- wavelength oscillations [16, 17, 1^ with Sm^^^ ^ 10^^^ eV^ and 



large mixing as a description of the solar neutrino deficit, it is instructive to determine how 
maximal or near-maximal solar and atmospheric vacuum oscillations may be accommodated 
within a three-neutrino universe. 

In this letter we use unitarity constraints to derive conditions on the three-neutrino mixing 
matrix under the assumption that both solar and atmospheric neutrinos undergo maximal 
mixing in vacuum. By maximal mixing we mean that the disappearance probabilities are 
equivalent to those for maximal two-neutrino mixing at the relevant mass scales. We find that 
there exists a unique three-neutrino mixing solution, up to trivial sign ambiguities, which 
conserves CP and has no oscillations of atmospheric Ue, consistent with indications from 
the Super-Kamiokande experiment 0. This mixing matrix predicts that solar Ue oscillate 
maximally into equal numbers of and Ur- We also investigate the consequences for three- 
neutrino mixing if the solar and atmospheric oscillations are not maximal. We find that for 
near-maximal solar and atmospheric neutrino mixing the — >■ Uf, and z/g — z/^ oscillation 
amplitudes for atmospheric and long-baseline experiments are approximately equal, and are 
limited by the deviation of the solar Ue disappearance amplitude from unity. Also, for any 
non-maximal solar and atmospheric neutrino mixing there must exist z/g Vt oscillations 
that may be visible in future long-baseline experiments. 

Oscillation probabilities. We begin our analysis with the survival probability for neutrino 
oscillations in a vacuum [|19| 

P{Ua -> Z/a) = 1 - 4 ^ PajPak siu^ Ajfc , (1) 

k<j 

where P^j = is the probability to find the a-fiavor state in the jth mass state (or vice 

versa), U is the neutrino mixing matrix (in the basis where the charged-lepton mass matrix 
is diagonal), Ajk = 6m],^L/AE = 1.27(5m2^/eV^)(L/km)/(E/GeV), ^m^^ = - ml, and 
the sum is over all j and k, subject to k < j. The matrix elements Uaj are the mixings 
between the fiavor (a = e, /x, r) and the mass (j = 1, 2, 3) eigenstates. We assume a neutrino 
mass spectrum with mass eigenvalues mi,m2 <C m^, where the solar oscillations are driven 
by A21 = Asun and the atmospheric oscillations are driven by A31 ~ A32 = A^tm- However, 
only the form of the neutrino mass matrix depends on the mass hierarchy assumption. All 
of our other conclusions require only the more general assumption that |As„„| = IA21I ^ 
IA31I ~ IA32I = |Aajm|, which also could apply, e.g., to the case < mi ~ m2, i.e., the 
solar neutrino vacuum oscillations are driven by a small mass splitting between the two larger 
masses. 
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The off-diagonal oscillation probabilites of this model are 



P{i^e = 4 P^sP^s sin^ A^tm - ARe{U,iU:^U;,U^2} sin^ A,„„ - 2 J sin 2A,,„ , (2) 

P(z/e ^Ur) = 4 P,sPr3 sm^ ^atm ' 4.Re{U eiUl^U^^U sin' A,,„ + 2 J sin 2A,,, , (3) 
PK^^^r) = 4P^3^'r3sin2A,t„-4i?e{f/^if/;2f/;if/r2}sin2A,,„-2Jsin2A,„,, (4) 

where the CP-violating "Jarlskog invariant" |2g is J = Jm {[/e2t/;3f^M2^M3 1 • The CP-odd 
term changes sign under reversal of the oscillating flavors. The solar amplitudes, quartic in 
the t/a/s, may themselves be expressed in terms of the Paj [0, |2^, as may the CP-odd 



amplitude J [|1|3 



We note that the CP- violating probability at the atmospheric scale is suppressed to 
order 5m^y„/5m^^^, the leading term having canceled in the sum over the two light-mass 
states. Thus, Piv^ vp) = P{i^/3 i^a) at the atmospheric scale. With MSW-enhanced 
solutions to the solar deficit, sin2Asun is effectively replaced by its average value of zero since 
^sun 3> 1 in this case. Only for the vacuum long-wavelength solution to the solar neutrino 
anomaly is there hope to measure the CP-odd term. 

We will use "amplitude" to denote the coefficients of the oscillating factors. The ampli- 
tudes of interest are: 

Alt = 4PeiPe2<(l-Pe3)', (5) 
A':L = 4P,3{P,l+P,2)=^P,3{l-P^3), (6) 
Alt = 4P,3(1-Pe3), (7) 
ACL = 4Pe3P^3<(l-Pr3)', (8) 
A'^L = 4P,,3P,3<(1-Pe3)', (9) 
All^ = 4P,3P.3<(1-Pm3)'- (10) 

In a three-neutrino model 

Aatm = A^^tm ~^ A'^tm ^ ^ud (11) 
^atm ~ ^atm ' ^atm ■ l-"-^/ 

The equality in Eq. @ results from the unitarity condition U*^Uf^2 + U*iUr2 = —U*^Ue2- 
The inequalities in Eqs. (H), (|8|)- (p!oD are obtained by maximizing the amplitude under the 
unitarity constraints 

Pe3 + P^3 + Pr3 = 1 , and (13) 
Pel + Pe2 + Pe3 = 1- (14) 

Thus, the amplitudes depend on just three independent P^/s, and the atmospheric ampli- 
tudes depend on just two of the probabilities chosen from Pe3, P^3, and Pt-s. 

Maximal solar and atmospheric oscillations. We first explore a solution for which both 
solar i/g and atmospheric data can be described by maximal oscillations at the relevant 
(5m^ scales, as suggested by the SuperK data P, ||. Maximal atmospheric oscillations require 
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Pfj.3 = ^- If solar oscillations are also required to be maximal, then Pes = and Pgi = Pt 



e2 



Up to trivial sign changes, unitarity then uniquely specifies the remaining elements of U : 
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(15) 



With the mixing matrix of Eq. 
are given by 



P{y„ 



(|15D , the off-diagonal oscillation probabilities in Eqs. (|D-(|D 



2 ^sun ; 

1 

2 



sin Ag^fi , 



(16) 
(17) 
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The zero t/es element implies that CP is conserved, J is zero, and the elements of U 
are real. It also implies that the only oscillations at the Aatm scale are z/^ Vr- The 
solar z/g oscillations are 50% into and 50% into Vri although the oscillation fiavors are 
not distinguishable because they have the same neutral current interactions and there is 
not enough energy to produce charge current events; however, detection of neutrinos from 
active galactic nuclei and gamma ray burst sources could distinguish the r-neutrino fiavor 
from Vf. and [|^. The mixing matrix for this bi-maximal scenario is similar to that 
obtained in Refs. and |2^, which give maximal solar oscillations {A^^^ = 1) and large 
but not maximal mixing of atmospheric neutrinos (A^fj^ = |). All of these scenarios also 
have Ue3 = 0, which implies maximal solar mixing and only oscillations at the /S.atm 

scale. 

A Majorana mass matrix or a hermitian Dirac mass matrix can be diagonalized by a single 
unitary matrix. For these two cases, we may therefore invert the process to obtain M in the 
flavor basis from the mass eigenvalues mj and the mixing matrix U . For the Majorana case 
Mmoj = UMdiagW^, while for the hermitian Dirac case MDirac = UM^iagU^ ■ The individual 

J2i iTT'jUajUg,, respectively. For real U, 



matrix elements are M„ 



— J2j ''^jUajUjSj and Maf3 — Z^j iiij<u aj<^ pji 

as here, there is no difference between these two cases. Reconstruction of the mass matrix 
m) gives 



from U in Eq. 



M = m 
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(19) 



in the (z/g, z/^. 



Ut-) basis, where the parameters of M are related to the mass eigenvalues as 

m = (20) 
e = (m2 + mi)/4m, (21) 
6 = ^2 (mi - m2)/4m . (22) 

Individual field redefinitions ipa —"ipa and ipj —ipj allow one to change the signs of 
individual rows and columns of M. We remind the reader that the form of M derived here 
applies in the basis where the charged-lepton mass matrix is diagonal. 
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The two relevant mass-squared differences are 



6ml^ ^^m\ 6ml^ = 8V26em\ (23) 

The atmospheric oscillation scale requires m ^ 0.035 eV. The vacuum long-wavelength solar 
solution then requires \6e\ ~ 10~®. Both 6 and e must be nonzero to generate a nonzero 
5m^„„. We reiterate that only the form of the mass matrix in Eq. (|TU]) depends on the 
assumption mi, 7712 <^ rris] the oscillation probabilites hold for the more general relation 

Near-maximal mixing. We next consider having both solar and atmospheric neutrino 
mixing close to maximal {A^f^ ^ 1 and A^f^^ ^ 1). Parametrizing the real mixing matrix in 
terms of Euler angles 

C1C3 — C1S3 Si 

S1C2C3 + S2S3 -S1C2S3 + S2C3 -C1C2 I I I , (24) 
-S1S2C3 + C2S3 S1S2S3 + C2C3 C1S2 
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Vr ) 








where Cj = cosOj and sj = sinOj, the maximal solar and atmospheric mixing of Eq. (15) 
corresponds to = 0, ^^2 = f ? and 6^ = j. Then for near-maximal mixing we take 

Ol=ei, 02 = ^ + £2, ^^3 = ^ + 63, (25) 

where the ej are small. The corresponding oscillation amplitudes are 

Alt ^ l-26?-4eL (26) 
A± ^ 1-Ael, (27) 

/t^e ^ ACT ^ 1 /le^ ~9p2 (HQ) 
^atm — ^atm — 2 a*™ — ^ ti • {^O) 

Here atmospheric disappearance is determined solely by €2- The solar z/g disappearance 
and atmospheric and long-baseline z/^ Ue and Ue oscillation amplitudes must satisfy 
the conditions 

2 aim, — -' '■atm — ■' '■atm ~ sun • \ J 

Hence for near-maximal mixing, z/^ — * Ue and — z^r oscillations in atmospheric and 
long-baseline experiments are small and approximately equal, and they are constrained by 
how close the solar mixing is to maximal. The atmospheric z/^ disappearance amplitude is 
independent of the other oscillation amplitudes in this limit. 

Non-maximal mixing. We also investigate the more general case of large but non-maximal 
solar neutrino mixing, which corresponds to nonzero Ue^, and large but non- maximal atmo- 



spheric neutrino mixing. A geometric expression for J pO| , equates its magnitude to 
the area of a triangle with sides which may be taken to be lUesU^^l, |f/e2^^t2|, and \UeiUf^i\. 
Since |f/e3| is known to be small, we use Eq.(|T5[) as an approximate guide for U and obtain 
l-'^l ~ |f^e3|/8 as the order of magnitude for the area of the triangle. This bound on \J\ is 
sufficiently small that we will set \J\ to zero and continue to work with real U and three 
independent parameters in the mixing matrix. If the oscillation amplitudes A^^^ and A'^f^ 
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are determined, e.g. 



Altm- Although pure z/„ ^ u, 



from experiment, then there remains only one free parameter, e.g., 
oscillations of atmospheric neutrinos are strongly disfavored, 
modest amounts of disappearance have not yet been ruled out. Therefore it is worthwhile 
to consider the effects of a nonzero Ues when the approximations of the previous section are 
no longer valid. 

Given the inputs v4^^„, v4^|^, and v4^f„, the three independent parameters of the matrix 
U can be determined up to some two-fold ambiguities, only one of which has phenomeno- 
logical consequences for vacuum oscillations. For example, knowing the atmospheric z/^ 
disappearance amplitude in Eq. determines U^^ = ±JPf^s: 



P, 



/x3 



It 



1 AfJ'l^ 
J- ^atm 



(30) 



where the — (+) means that u^^ is more closely associated with z/2 (z/3). Knowing the atmo- 
spheric z/g disappearance amplitude in Eq. (0) determines Ues = Ty/Pes'. 



e3 



1-A 



atm 



(31) 



In Eq. (^T]) we have chosen the solution which gives Pes < \ since Pe3 > \ implies that the 
solar z/g disappearance amplitude must be less than \, which for the long- wavelength 



solar solution is disfavored by the data. Then the unitarity relation in Eq. (|T3D yields 



1 _ /le^ _|_ ./1 _ A^^l^ 



^atm 



where the ± in Eq. ( p2D is correlated with the =F in Eq. 
relation in Eq. (|I^) imply 



(32) 



Also, Eq. (H) and the unitarity 



Pel 
Pe2 



l-Pe3 + V(l-^e3)'-A 
I -Pes- 



eq 
sun 



(1 — Pe3)^ ~ Astn 



(33) 
(34) 



where we have assumed that is more closely associated with z/i than with z/2 (the other 
possibility, that z/g is more closely associated with z/2, has similar vacuum oscillation phe- 
nomenology). In principle, there is a constant term in the solar oscillation probability 
2Pe3(l — Pes) = |/lafm) f^om averaging of the sin"^ A atm term. However, if A^f^ is small, as 
suggested by the non-observation of z/g oscillations in the atmospheric data, this term can be 
ignored to first approximation The remaining elements of U may also be determined 
from unitarity (e.g., by using the values of P^s, Pes, Prs, Pei, and Pe2 and the form of Eq. ( p4|) 
to determine the three Euler mixing angles). 

Once the oscillation parameters are fixed (up to the two- fold ambiguity in Eq. (|30D ), 
other oscillation amplitudes may be calculated; for example, Eqs. ( pUD and (^TJ) combined 
with Eq. (I) yield 

'1 - yJl-Alt] [l T Vl-^ail , (35) 



.fie 
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where the =F is correlated with the =F in Eq. (|30|). Contours of A'^^^ in the {^A^jlrni-^atm) plane 



are shown in Figs. 1(a) and 1(b). Once J^f^rn and J^atm are measured, the value of A^^^ is 
predicted. 

In principle the three inputs can be determined by the solar and atmospheric neutrino 
data. For instance, A^^^ is directly determined from vacuum long-wavelength fits to the 
solar neutrino spectrum. Also, the two independent atmospheric oscillation amplitudes can 
be found from measurements of 



N.INI = {\-<S> Alt)^T <S> A^J, 



atm y 



(36) 



and 



'1-<S> A 



atm. 



+ r ^ < S > A^^^ , 



where N° and A*"" are the expected number of atmospheric e and /i events, r = N°/N° 
and < 5* > is sin^ A^tm averaged over the appropriate atmospheric neutrino spectrum. 



(37) 



M 2' 



The existence of a real solution to Eqs. 



and 



requires 



Pe3<l-\/A 



(38) 



The limit in Eq. (pq) becomes stringent when is close to unity. Current data suggest 
Agun — 0-6; which gives Pes < 0.23 and hence A^f^ < 0.7. These limits will improve as more 
solar and atmospheric data become available. As a check, note that in Eq. ( |5BD maximal 
mixing of solar neutrinos (A^fn = 1) automatically leads to jf/esl = ^aim = ^alm = 0; i-^-, 
all atmospheric oscillations are z/^ — > z/^. From our approximate bound for \J\ when Ues is 
small, we obtain \J\ ^ 0.06, which verifies that CP- violation effects will be small. 

Another limit on P^^ comes from the CHOOZ reactor experiment [^] that measures Ue 
disappearance 

(39) 

I) gives the 



/.e^ < n 9 



which applies for Sml^^ ^ 2 x 10 
bound 



eV^. When combined with Eq. (|3lD, Eq 
Pes < 0.05 . 



(40) 



However, for Sml^^ < 10 



eV^ there is no limit at all. The proposed experiment to detect 



reactor anti-neutrinos in the BOREXINO detector ||29| can test a value of A'^^^ as low as 0.4 
(corresponding to a limit Pes < 0.18) for 5^^^^ ~ 10"^ eV^. 

The new physics predicted by the non-maximal mixing model is Ue Vt oscillations with 
leading probability (i.e., ignoring oscillations involving Asnn) 



4 PesPrS Sin^ /\atm , 



(41) 



which from Eqs. (0), (11), and ([To|) has an oscillation amplitude 



/ler 

atm 



1-A 



atm 



Aatm =t A/ 1 



■'^atm 



(42) 



where the ± correlates with the =F in Eq. (^ 

once and A^f^ are known, modulo the two-fold ambiguity in Eq. 



As with A^l^, ^alm is completely determined 



Contours of ^^[m. 
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in the {A'^tm^^ltm) plane are shown in Figs. 1(c) and 1(d). We note that the approximate 
equahties involving the atmospheric oscillation amplitudes in Eq. (|29|) can also be derived 
by taking the Af^t, 1, Alf^ limit in Eqs. (|3|) and (||). 

The z/e — »• Ut oscillations could be observed by long-baseline neutrino experiments with 
proposed high intensity muon sources pO, [TO, [1^, which can also make precise measurements 



of z/g, i^^ i^T, and z/g — * oscillations. Sensitivity to yl^[„(5m^j„/eV^)^ > 2.5 x 

-9 



10~ is expected for the parameter ranges of interest here. The combination of all 
measurements (solar, atmospheric, and long-baseline) would serve as a consistency check on 
the model. 



Other consequences. Because the highest mass scale is J5m^i^ ^ 0.07 eV, there is 



little neutrino contribution to hot dark matter in this model ||31| , p2| . Such a mass can be 
generated for Majorana neutrinos by the see-saw mechanism with a heavy Majorana scale 
of about 5 X 10^^ GeV (radiative corrections may reduce this value by about a factor of two 
||33|| ) if the Dirac mass is of order the top mass. Also, if mi <C m2 < in this model, then 
1712/1713 ~ 10~^ is approximately equal to {nic/mt)'^, as predicted in some see-saw models 
| |33| |. The small value of the mass matrix element Mi,^,y^ ^ 10^^ eV means that the model 
has an unmeasurably small contribution to neutrinoless double-beta decay p3. The small 



values of t/es and mean that there are no observable consequences for tritium beta-decay 



35| . However, the heaviest mass does provide a relic neutrino target for a mechanism that 
may generate the cosmic ray air showers observed above ^ 10^° eV . 

Discussion. More precise measurements of the spectrum shape and day-night effects 
by SuperK and the Sudbury Neutrino Observatory (SNO) |^ could determine if the solar 
oscillations are vacuum long-wavelength. Measurements of a seasonal variation of the flux 
of solar ^Be and pep neutrinos , which can be measured in the BOREXINO experiment 
p9| , would confirm that the solar oscillations are vacuum long-wavelength. The HELLAZ 
experiment plans to measure the solar pp neutrino spectrum, which would also be useful 
in establishing the vacuum long- wavelength solution. All of these solar experiments should 
also be able to pin down the value A^f^^ in the vacuum long- wavelength scenario. 

Further measurements of atmospheric neutrinos will more precisely determine the am- 
plitudes A'^f^ and A^f^, completely fixing (modulo a two-fold ambiguity) all oscillation 
amplitudes in atmospheric and long-baseline experiments. The MINOS experiment is 
expected to test for z/^ — ^ z/g and z/^ — > Ur oscillations for Sm^^^ > 10"^ eV^. Together 
the solar and atmospheric measurements could put strong limits on Ug^, which governs the 
amount of z/g — >■ mixing that could be seen in future long-baseline experiments such as 



those discussed for the proposed muon collider at Fermilab [Q. For near-maximal solar and 



atmospheric neutrino mixing the z/^ — > and z/g i^t oscillation amplitudes for atmospheric 
and long-baseline experiments are approximately equal, and are limited by the deviation of 
the solar z/g disappearance amplitude from unity. 
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Figure 1: Correlations of atmospheric oscillation amplitudes: shown are contours of Aatm 
the {A'^f^,A'^tm) plane for the (a) negative root, and (b) positive root solutions for P^^ in 
Eq. (0), which correspond to being more closely associated with 1^2 and z/3, respectively. 
Also shown are contours of A'^^ for the (c) negative root, and (b) positive root solutions. 
The point with bi-maximal solar and atmospheric neutrino oscillations is noted. 
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